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Nuclear reactions generate energy in nuclear reactors, in stars, and are responsible for the ex- 
istence of all elements heavier than hydrogen in the universe. Nuclear reactions denote reactions 
between nuclei, and between nuclei and other fundamental particles, such as electrons and photons. 
A short description of the conservation laws and the definition of basic physical quantities is pre- 
sented, followed by a more detailed account of specific cases: (a) formation and decay of compound 
nuclei; (b)direct reactions; (c) photon and electron scattering; (d) heavy ion collisions; (e) forma- 
tion of a quark-gluon plasma; (f) thermonuclear reactions; (g) and reactions with radioactive beams. 
Whenever necessary, basic equations are introduced to help understand general properties of these 
reactions. Published in Wiley Encyclopedia of Physics, ISBN-13: 978-3-527-40691-3 - Wiley- VCH, 
Berlin, 2009. 
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I. INTRODUCTION 

The collision of two nuclei can give place to a nuclear 
reaction and, similarly to a chemical reaction, the final 
products can be different from the initial ones. This pro- 
cess happens when a target is bombarded by particles 
coming from an accelerator or from a radioactive sub- 
stance. It was in the latter way that Rutherford observed, 
in 1919, the first nuclear reaction produced in laboratory, 

a+ V'N^ 'lO +p, (1) 

using a-particles from a ^^^Bi sample. 

As in eq. [TJ other reactions were induced using a- 
particles, the only projectile available initially. With the 
development of accelerators around 1930, the possibili- 
ties multiplied by changing the energy and mass of the 
projectile. Today it is possible to bombard a target with 
protons of energy greater than 1 TeV (1 TeV = 10^^ eV 
= 1.602x 10^^ joules), and beams of particles as heavy as 
uranium are available for study of reactions with heavy 
ions. 

Sometimes we can have more than two final products 
in a reaction, as in 

p + ^''N^ ^Be + 2a, 
p + ^^Na^ 2^Ne +p + n, (2) 

or just one, as in the capture reaction p -I- ^''Al — ^®Si*, 
where the asterisk indicates an excited state, which usu- 
ally decays emitting 7-radiation. Under special circum- 
stances, more than two reactants is possible. Thus, for 
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example, the reaction a+ a+ a — )■ ""^^C can take place 
in the overheated plasma of stellar interiors. 

The initial and final products can also be identical. 
This case characterizes a process which can be elastic, as 
inp + ^^O — p + ^^O, where there is only transfer of 
kinetic energy between projectile and target, or inelastic, 
as in n + ^^O — )■ n + ^^O*, where part of the kinetic 
energy of the system is used in the excitation of ^^O. 

Naturally, nuclear reactions are not limited to nuclei. 
They could involve any type of particle, and also radia- 
tion. Thus, the reactions 

7 + ^^Cu ^ ^^Ni + p, 

^ + 233u ^ 90Rb + "ics + 2n, (3) 

are examples of nuclear processes induced by gamma ra- 
diation. In the first case a 7-ray knocks a proton off 
^■^Cu and in the second it induces nuclear fission of ^'^■^U, 
with the production of two fragments and emission of two 
neutrons. 

Unlike a chemical reaction, the resulting products of a 
nuclear reaction are not determined univocally: starting 
from two or more reactants there can exist dozens of final 
products with an unlimited number of available quantum 
states. As an example, the collision of a deuteron with 
^■^^U can give place, among others, to the following reac- 
tions: 

d + 238u ^ 240Np + 
-> 239p^p _^ 
239U + 

^ 237u + t. (4) 

In the first of them the deuteron is absorbed by the ura- 
nium, forming an excited nucleus of ^'^''Np that de-excites 
by emitting a 7-ray. The two following are examples of 
stripping reactions, in which a nucleon is transferred from 
the projectile to the target. The last one exemplifies the 
inverse process: the deuteron captures a neutron from 
the target and emerges out as (tritium). This is de- 
noted as a pick-up reaction. Another possibility would 
be, in the first reaction, that ^'^'^Np fissions instead of 
emitting a 7-ray, contributing with dozens of possible fi- 
nal products for the reaction. 

Each reaction branch, with well defined quantum 
states of the participants, is known as channel. In [4j 
for the entrance channel d + ^"^^U, there are four pos- 
sible exit channels. Notice that a different exit channel 
would be reached if some of the final products were in 
an excited state. The probability that a nuclear reaction 
takes place through a certain exit channel depends on 
the energy of the incident particle and is measured by 
the cross section for that channel. 

Nuclear reactions proceed through many possible dis- 
tinct mechanisms. For instance, in direct reactions the 
projectile and the target have an interaction of short du- 
ration, with possible exchange of energy or particles be- 
tween them. Another mechanism involves the fusion of 
the projectile with the target, the available energy being 



distributed among all nucleons, forming a highly excited 
compound nucleus. The decay of the compound nucleus 
leads to the final products of the reaction. 

In high energies collisions, the nuclei fragment and 
particles that were not initially present are produced 
(for instance, pions, kaons, etc.). The reactions proceed 
through an intermediate phase in which the nuclear mat- 
ter is compressed. At very high energies the quarks and 
gluons inside nucleons become "free" for a short time 
forming the quark-gluon plasma. The study of high en- 
ergy reactions with nuclei is very important for a better 
understanding of what happens during spectacular stel- 
lar explosions (supernovae) and in the interior of compact 
stars, as for instance, neutron stars. The study of nuclear 
reactions at high energies allows to obtain information on 
the equation of state of nuclear matter. 

II. BASIC PRINCIPLES 
A. Conservation laws 

Several conservation laws contribute to restrict the 
possible processes in a nuclear reaction. 

I- Baryonic number - There is no experimental ev- 
idence of processes in which nucleons are created or 
destroyed without the creation or destruction of corre- 
sponding antinucleons. The application of this principle 
to low energy reactions is still more restrictive. Below the 
threshold for the production of mesons (~ 140 MeV), no 
process related to the nuclear forces is capable to trans- 
form a proton into a neutron and vice-versa, and pro- 
cesses governed by the weak force (responsible for the /3- 
emission of nuclei) are very slow in relation to the times 
involved in nuclear reactions (~ 10~^^ to 10~"'^^s). In this 
way, we can speak separately of proton and neutron con- 
servation, which should show up with same amounts in 
both sides of a nuclear reaction. 

2 - Charge - This is a general conservation principle 
in physics, valid in any circumstance. In purely nuclear 
reactions it is computed making the sum of the atomic 
numbers, which should be identical, at both sides of the 
reaction. 

3 - Energy and linear momentum - These are two of the 
most applied principles in the study of the kinematics 
of reactions. Through them, angles and velocities are 
related to the initial parameters of the problem. 

4 - Total angular momentum - is always a constant of 
motion. In the reaction 

i"B + %e^ + "C, (5) 

^"B has / = 3 in the ground state, whereas the a-particle 
has zero angular momentum. If it is captured in an s- 
wave {k = 0), the intermediate compound nucleus is in 
a state with Ic = 3. Both final products have intrinsic 
angular momenta equal to 1/2. Hence, their sum is or 
1. Therefore the relative angular momentum of the final 
products will be = 2, 3 or 4. 
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5 - Parity - is always conserved in reactions governed 
by the nuclear interaction. In the previous example, ^''B, 
■^He and the proton have equal parities, while ^'^C has odd 
parity. Therefore, if h = 0, we necessarily have If — 3. 
Thus, the orbital momentum of the final products of eq. 
[5] is determined by the joint conservation of the total 
angular momentum and of the parity. 

6 - Isospin - This is an approximate conservation law 
that is applied to light nuclei, where the effect of the 
Coulomb force is small. A nuclear reaction involving 
these nuclei not only conserves the z-component of the 
isospin (a consequence of charge and baryonic number 
conservation) but also the total isospin T. Reactions that 
populate excited states not conserving the value of T are 
strongly inhibited. An example is the reaction d + ^^O 

a+ "N, where the excited state 0+, with 2.31 MeV, 
of is about a hundred times more populated than 
the 1^ ground state. Conservation of energy, angular 
momentum and parity do not impose any prohibition for 
that channel, whose low occurrence can only be justified 
by isospin conservation: the ground states of the four 
participant nuclei in the reaction have all T = and the 
state 1+ of "N has T = 1. 




FIG. 1: (Top) Rutherford's scattering experiment with a- 
particles scattering on a gold foil. (Bottom) Microscopic in- 
terpretation of the experimental results. 



B. Kinematics 

We consider a typical reaction, in which the projectile 
a and the target A gives place to two products, b and 
B. This can also be expressed in the notation that we 
used so far, a + A — > b + B, or even in a more compact 
notation, A(a,b)B. 

Often, a and b are light nuclei and A and B, heavy 
ones; the nucleus b being emitted at an angle 9 and its 
energy registered in the laboratory system. The recoiling 
nucleus B has usually a short range and cannot leave the 
target. It is convenient to introduce the Q-value of the 
reaction which measures the energy gained (or lost) due 
to the difference between the initial and final masses: 

Q ^ {rria + m,A - rub - mB)c'^ ■ (6) 

Using energy and momentum conservation in the reac- 
tion, one gets: 

\ TUB J V ™B/ 

2 

- \/ niambEaEb cos 9. (7) 

rriB 

From this relation one obtains that, when Q is nega- 
tive, an energy threshold exists for the incident particle, 
Et, as a function of the angle 0, below which nuclei b are 
not observed in that angle: 

^ _ -QniBirriB + rnt) 

maTOftCOS^fl + {niB + mi,){mB — rUa) 

This example shows the power of conservation laws in 
the analysis of nuclear reactions. Similar to this example. 



very useful relations can be derived in the relativistic 
regime, and/or using other conserved quantities. 

C. Cross sections 

Figure [T] shows schematically a typical scattering ex- 
periment. In fact this is the sketch of E. Rutherford 
experiment in 1910 [5]. Projectiles (here, a-particles) 
from a source pass through a collimator and collide 
with a target (gold foil). Some projectiles are scattered 
by the target and reach the detector (here, a fluores- 
cent screen). Rutherford expected them to go straight 
through the target, with perhaps minor deflections. Most 
did go straight through, but to his surprise some particles 
bounced straight back! Rutherford hypothesized that the 
positive alpha particles had hit a heavy mass of positive 
particles, which he called the nucleus. 

Let us consider an experiment as that represented in 
figure [2] measuring the count rate of events leading to the 
population of channel- a, Na{^l, Assuming that in- 

teractions between beam particles can be neglected, the 
count rate A^Q(r2, AfJ) should be proportional to the in- 
cident flux J, to the fraction of solid angle Ail where the 
particles are scattered to, and to the number of target 
particles per unit volume. Hence, Na{^^ Ail) can then be 
written as Na{^,/\n) = {An ■ n ■ J) daa{^)/dVl. The 
constant of proportionality, 

rfao(») ^ iVa(»,A») 

dn An-n-j ' ^ ' 

is called the differential cross section for channel-a. This 
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FIG. 2: Schematic representation of a scattering experiment. 



quantity is very useful, since it does not depend on exper- 
imental details (detector size, incident flux, target thick- 
ness) . It depends exclusively on the physics of the pro- 
jectile and target particles. 

Frequently, one is interested in the angle integrated 
total cross section given by 



1 dn 






dfl 



(10) 



In nuclear physics cross sections are usually measured 
in units of barns. 1 barn corresponds to an area of a circle 
of diameter approximately equal to 8 fm = 8 xlO~^^ m. 
For the comparison between theory and experiment, it 
is necessary to have cross sections in the same reference 
frame. Of course, the measured cross section is obtained 
in the laboratory frame, where the target is at rest. From 
the theoretical point of view, however, it is important 
to take advantage of the translational invariance of the 
projectile-target Hamiltonian, introducing the center of 
mass frame (CM). The transformation is done using the 
laws of conservation of energy and momentum. 

The Rutherford experiment has a simple classical de- 
scription in which one assumes that the a-particles follow 
hyperbolic trajectories in the Coulomb field of the heavy 
target. The scattering at an angle 9 can be easily calcu- 
lated and with that one is able to calculate the Coulomb, 
or Rutherford cross section. The result is 
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i^^J '''''' I2J ' 



(11) 



where Zp (Zt) is the projectile (target) charge number 
and E is the CM energy. 



D. Elastic scattering 

When a beam of particles, represented in quantum me- 
chanics by a plane wave, hits a nucleus the wave function 
is modified by the presence of a scattering potential V{r), 
responsible for the appearance of a phase in the outgo- 
ing part of the wave. Elastic scattering is just one of the 
channels for the which the reaction can proceed and is 



known as elastic channel. Inelastic scattering and all the 
other channels are grouped in the reaction channel. 

The occurrence of a nuclear reaction through a given 
reaction channel leads to a modification of the outgoing 
part of the wave function not only by a phase factor, but 
also by changing its magnitude, indicating that there is 
a loss of particles in the elastic channel. For a projectile 
with momentum p = hk, this can be expressed by 
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1 

2i ^ 

1=0 



{21 + Pi{cos6 



i(kr-liT/2) -i{kr~l-!i/2) 



kr 



(12) 

where the complex coefficient rji is the factor mentioned 
above. If 77 = 1, the sum in eq. [l2]can be done analyt- 
ically, leading to ~ exp(ik.r), i.e. a plane wave. But 
if rji — exp[i(5;], with Si real, the incoming and outgoing 
waves have the same magnitude, i.e. the scattering is 
elastic. 



The sum in eq. 12 is know as partial wave expansion of 
the scattering wave. The label Z = 0, 1, 2, • • • = (s, p, d, 
• • • waves) denotes the contribution of a particular angu- 
lar momentum (in units of h) to the total wavefunction. 
Classically the angular momentum of an incident particle 
is given by I = kR, where R is known as impact param- 
eter which is the perpendicular distance to the target if 
the projectile were undeflected. In quantum mechanics I 
is not continuous, varying in steps of one. 

We can rewrite eq. [12] as a sum of a plane wave and 
a scattering outgoing wave in the form ^ ~ exp(ik.r) -|- 
f{9)e'^^'^/r, where f{0) accounts for the distortion of the 
outgoing wave at the scattering angle 9. It is known as 
the scattering amplitude: 



/(^) = ^ E(2^ + - l)fl(cos0)- (13) 



(=0 



One can now calculate the cross section by counting the 
number of particles that are scattered through angle 9. 
This can be done by calculating the particle current asso- 
ciated with the wavefunction One gets the differential 
scattering cross section 
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^(2? + l)(l-77i)Pi(cos6 



The total scattering cross section, eq. [TO] becomes 



cr^ =7rA2^(2/ + l)|l- 



(14) 



(15) 



with A = A/27r = 1/fc. 

E. Reaction cross sections 

To calculate the reaction cross section it is necessary 
to compute the number of particles that disappear from 
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FIG. 3: Elastic scattering angular distribution of protons of 
30 MeV on Ca, Sn and Pb. The curves are adjusted to the 
experimental data of ref. [3] 



For Tji = 0, the scattering and reaction cross sections are 
identical, yielding the total cross section 

cr = (jr + (Je 7r(i? + A)^ + 7r(i? + A)^ = 27r(i? + A)^ (18) 

that is twice the geometric cross section! 

The presence of the scattering part, that turns the re- 
sult [is] apparently strange, can be interpreted as the ef- 
fect of diffraction of the plane waves at the nuclear sur- 
face. This effect leads to "shadow" behind the nucleus 
decreasing its apparent diameter so that, at a certain dis- 
tance, the perturbation caused by the presence of the nu- 
cleus disappears and the plane wave is reconstructed. In 
this situation we can say that the part of the beam which 
is diffracted has to be the same as the part that is ab- 
sorbed, justifying the equality of ar and Ue. The diffrac- 
tion phenomenon appears clearly in the elastic scattering 
or inelastic angular distribution (differential cross section 
as function of the scattering angle): figure |3] exhibits an- 
gular distributions for the elastic scattering of 30 MeV 
protons on *°Ca, ^^°Sn and ^°^Pb. The oscillations in the 
cross sections are characteristic of a Fraunhofer diffrac- 
tion figure, similar to light scattering by an opaque disk. 
The angular distance A9 between the diffraction minima 
follows closely the expression A9 = h/pR, characteristic 
of diffraction phenomena. 



the elastic channel, what is measured by the flux of the 
current of probability vector through a spherical surface 
of large radius centered at the target, calculated with the 
total wave function of eq. 



12 One finds 



= ^P^(2; + i)(i-b|2). 

1=0 



(16) 



From equations 15 and 16 when [Ty/j = 1 the reaction 
cross section is zero and we have pure scattering. The 
contrary, however, cannot happen, as the vanishing of ae 
also implies in the vanishing of ar- In general there is 
a region of allowed values of rji for which the two cross 
sections can coexist. 

The maximum of tr^ happens for rji ^ 0, what cor- 
responds to total absorption. Let us suppose that the 
absorption potential is limited to the surface of a nucleus 
with radius i? S> A, that is, that all the particles with im- 
pact parameter smaller than the radius R are absorbed. 
That is equivalent to say that all particles are absorbed 
for Z < i?/A. In this case 



J?./ A 

7rA2^(2? + l) = 7r(i? + A)^ 



(17) 



This is the value that would be intuitively adequate for 
the total cross section, i.e., equal to the geometric cross 
section (the part A can be understood as an uncertainty 
in the position of the incident particle). But, we saw 
above that the presence of scattering is always obligatory. 



F. Excitation functions 

When the projectile has a very low energy, fc — 0, 
and in particular, I — kR ^1. As an example con- 
sider the scattering of neutrons with I = and ignore the 
spins of the neutron and of the target. In this case the 
Schrddinger equation for the radial motion of the neutron 
with respect to the target is 



dr^ 



+ k^uo = (r > R). 



(19) 



This is valid for the radial wave function uq at distances 
r larger than the channel radius R = Ra + Ra, with Ra 
and Ra being the radii of the projectile and of the target, 
respectively. The solution of eq. [19] is 



— ikr 



(r > R). 



(20) 



A radial wave function inside the nucleus should connect 
to the external function [20] with a continuous function 
and its derivative at r = _R. That is, the function 



fi = R 



dui I dr 



ui 



(21) 



r-i?. 



must have identical values if calculated with the internal 
or the external function and this condition creates a re- 
lationship between /; and ry/. Hence, the knowledge of /; 
leads to the knowledge of the cross sections. For neutrons 



with / = 0, the application of 20 results in 



/o = ikR 



m 



-2ikR 



VO 



-2ikR ■ 



(22) 
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from which we extract 



^ .fo±ikR^^2ikR 
fo - ikR 



(23) 



If /o is a real number, then |7yoP = 1- The reaction 
cross section, eq. |16[ will be zero and we have pure scat- 
tering. 

Using eq. [23) the scattering cross section, eq. [15] can 
be written as 



Apot I ^ J 



(24) 



with Ar 



exp(2ifci?) — 1. The sepa- 
ration of the cross section in two parts has physical justi- 
fication: Apot does not contain /o and therefore does not 
depend on conditions inside the nucleus. It represents 
the situation where the projectile does not penetrate the 
nucleus as in the idealized situation where the nucleus is 
considered an impenetrable hard sphere. The wave func- 
tion is zero inside the nucleus and uq vanishes at r = _R, 
implying /g — >■ oo and A^es 0. Hence, Ap^t is the only 
responsible for the scattering. 



Inserting eq. 23 in [16^ and using 



fo = fR + -if I 



we have 



-AkRfi 



/I + (// - kRY 



(25) 



(26) 



an equation that is useful when we study the presence of 
resonances in the excitation function (cross section as a 
function of the energy) . 



III. STATISTICAL REACTIONS 

A. Compound nucleus 

When a low energy neutron (< 50 MeV) enters the 
range of nuclear forces it can be scattered or begin a se- 
ries of collisions with the nucleons. The products of these 
collisions, including the incident particle, will continue in 
their course, leading to new collisions and new changes 
of energy. During this process one or more particles can 
be emitted and they form with the residual nucleus the 
products of a reaction that is known as pre- equilibrium. 
At low energies, the largest probability is the continua- 
tion of the process so that the initial energy is distributed 
among all nucleons, with no emitted particle. The final 
nucleus with A + I nucleons has an excitation energy 
equal to the kinetic energy of the incident neutron plus 
the binding energy the neutron has in the new, highly un- 
stable, nucleus [4j. It can, among other processes, emit 
a neutron with the same or smaller energy to the one 
absorbed. The de-excitation process is not necessarily 
immediate and the excited nucleus can live a relatively 
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FIG. 4; Cross sections for tlie reactions shown in 1271 The 
scales of the upper axis (energy of the protons) and lower axis 
(energy of the a-particle) were adjusted to correspond to the 
same excitation energy of the compound nucleus [5] . 



long time. We say that there is, in this situation, the 
formation of a compound nucleus as intermediary stage 
of the reaction. In the final stage the compound nucleus 
can evaporate one or more particles, fission, etc. In our 
notation, for the most common situation in which two 
final products are formed (the evaporated particle plus 
the residual nucleus or two fission fragments, etc.) we 
write: 



A-^ C* 



B 



the asterisk indicating that the compound nucleus C is 
in an excited state. 

The compound nucleus lives long enough to "forget" 
how it was formed and the de-excitation to the final prod- 
ucts b and B only depends on the energy, angular momen- 
tum and parity of the quantum state of the compound 
nucleus. An interesting experimental verification was ac- 
complished by S. N. Ghoshal in 1950 [5]. He studied 
two reactions that take to the same compound nucleus, 
^^Zn*, and measured the cross sections of three different 
forms of decay, as shown below: 



^^Cu 

6"Ni 



/ 632n +n 

\ 622n + 2n (27) 



If the idea of the compound nucleus is valid and if one 
chooses the energy of the proton and of the incident 
a-particle to produce the same excitation energy, then 
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the cross section for each one of the three exit channels 
should be independent of the way the compound nucleus 
is formed. That is, the properties of the compound nu- 
cleus do not have any relationship with the nuclei that 
formed it. This is confirmed in figure [4j where one sees 
clearly that the cross sections depend practically only on 
the exit channels. 

The angular distribution of fragments, or evaporated 
particles, of a compound nucleus should be isotropic in 
the center of mass, and this is verified experimentally. 
However, the total angular momentum is conserved and 
cannot be "forgotten". Reactions with large transfer of 
angular momentum, as when heavy ions are used as pro- 
jectiles, can show a non-isotropic angular distribution in 
the center of mass system. 

The occurrence of a nuclear reaction in two stages al- 
lows the cross section for a reaction A(a,b)B to be written 
as the product, fT(a,b) = crc(a,A)P(b), where ac{a,A) is 
the cross section of formation of the compound nucleus 
starting from the projectile a and the target A and P(b) 
is the probability that the compound nucleus emits a 
particle b leaving a residual nucleus B. If the quantum 
numbers of entrance and exit channels are well specified, 
i.e., if the reaction begins at an entrance channel a and 
ends at an exit channel /3, one can write 



a(a,/3) =a,(a)P(/3). 



(28) 



We can associate the probability P{l3) to the width 
of the channel /3 and write: 



(29) 



where F is the total width, that is, r = h/T is the half- 
life of disintegration of the compound nucleus. Eq. |29] 
just expresses the fact that the decay probability through 
channel /3 is the decay rate through that channel divided 
by the total decay rate. In the competition between the 
several channels /3, the nucleons have clear preference 
over the 7-radiation whenever there is available energy 
for their emission and among the nucleons the neutrons 
have preference as they do not have the Coulomb barrier 
as an obstacle. Thus, in a reaction where there is no 
restriction for neutron emission we can say that 



F. 



(30) 



where F„ includes the width for the emission of one or 
more neutrons. 

The study of the function P(/3) is done in an evapora- 
tion model that leads to results in many aspects similar to 
the evaporation of molecules of a liquid, with the energy 
of the emitted neutrons having the form of a Maxwell- 
Boltzmann distribution 



I{E) cx Eexp(-j^dE 



(31) 



with / measuring the amount of neutrons emitted with 
energy between E and E + dE. The quantity 6, with di- 
mension of energy, has the role of a nuclear temperature. 




FIG. 5: Energy spectrum of neutrons evaporated by a com- 
pound nucleus. 



It is related to the density of levels oj of the daughter 
nucleus B by 



with 



1 _ dS 
9~dE' 



5 = lnw(i;) 



(32) 



(33) 



where dS/dE is calculated for the daughter nucleus B at 
the maximum excitation energy that it can have after the 
emission of a neutron. That is, in the limit of emission 
of a neutron with zero kinetic energy. 

The level density uj{E) is a measure of the number of 
available energy states for the decay of the compound 
nucleus in the interval dE around the energy E. In this 
sense, the relationship [33] is, neglecting the absence of 
the Boltzmann constant, identical to the thermodynamic 
relationship between the entropy S and the number of 
states available for the transformation of a system. Eq. 
|32| is the well-known relation between the entropy and 
the temperature. 

A simple expression for the energy dependence of 
the state density is provided by the equidistant spac- 
ing model which assumes that the one-particle states are 
equally spaced with spacing d, and that the total energy 
of the nucleus is simply obtained by adding the energies 
of the constituent nucleons. The solution of this problem 
can be obtained from statistical mechanics [lOj : 

p [E) ^ cxp I 

Extensive analyzes of experimental data show that for 
nuclei far from the region of the magic nuclei a varies 
linearly with A (or with N and Z), as shown in fig. [6j 



A 



MeV" 



(34) 



It is found that k ^ 7.5 
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FIG. 6: Values of the level density parameter a as a func- 
tion of the neutron number. Far from magic regions, a « 
A^/4.5 MeV~^ which approximately corresponds to o = A/7.5 
MeV-i [ig. 
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FIG. 7: Total cross section for low energy neutrons hitting 
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B. Energy spectrum of neutrons 



Let us now assume that there is only elastic scattering 
or formation of compound nucleus. We further assume 
that the elastic scattering is purely potential, without res- 
onant elastic scattering. That is, there is no re-emission 
of neutrons with energy equal to the incident projectile. 
This is equivalent to say that the probability that the 
exit channel is the same as the entrance channel is very 
low. According to these assumptions, the cross section 
of eq. 28 is the reaction cross section of eq. 16 



We can still write the wave function inside the nucleus 
as just an incoming wave 



{r < R), 



(35) 



where K = ^2m(E~ Vo)/fi. is the wavenumber inside 
the nucleus, and it is assumed that the neutron with total 
energy E is subject to a negative potential Vq- Eq. [35] is 
clearly a crude simplification when the incident neutron 
interacts in a complicated way with the other nucleons 
in the nucleus. It allows, however, to explain the average 
behavior of the cross sections for low energies. Starting 



with eqs. 21 and 35 one determines the value of /q: 



/o = -iKR 



and from eq. [26) we get 



{k + Kf 



(36) 



(37) 



for the cross section of compound nucleus formation for 
neutrons with I — Q. At low energies, E « |Vb|, thus 
k « K. Under these conditions, CTc = 4:TT/kK, since 
A = k^^ » R. Thus, (Tc varies with 1/k. That is, 



1 



CTn DC 



(38) 



where v is the velocity of the incident neutron. This is 
the well-known 1/v law that governs the behavior of the 
capture cross section of low energy neutrons. Figure [7] 
exhibits the excitation function (cross section as function 
of the energy) for the reaction n-l-^'^^U. The cross section 
decays with 1/v up to 0.3 eV, where a series of resonances 
appear. 



The energy distribution of neutrons emitted by a com- 
pound nucleus has the aspect of the curve shown in figure 



[5] Only the low energy part obeys 31 and the reason is 
simple: the emission of a low energy neutron leaves the 
residual nucleus with a large excitation energy, and the 
level density is very high. The large density of final states 
turns the problem tractable with the statistical model 
that leads tojSl] In the opposite situation are the low en- 
ergy states of the residual nucleus. These isolated states 
appear as peaks in the tail of the distribution. When the 
emission is of a proton or of another charged particle, the 
form of figure [3l] is distorted, the part of low energy of 
the spectrum being suppressed partially by the Coulomb 
barrier. 



C. Resonances 

To understand why there are resonances, we shall use 
again the simple model of a single particle subject to a 
square-well potential. We know that inside the well the 
Schrodinger equation only admits solution for a discrete 
group of values of energy, Ei, -E2 , • ■ • ■ A particle is con- 
fined to the interior of the well by reflections that it has 
at the surface of the well. In these reflections the wave 
that represents the particle should be in phase before 
and after the reflection and this only happens for a finite 
group of energies. Outside the well the Schrodinger equa- 
tion does not impose restrictions and the energy can have 
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any value. But we know, from the study of the passage 
of a beam of particles through a potential step, that the 
discontinuity of the potential at the step provokes reflec- 
tion even when the total energy of the particles is larger 
than the step, a situation where classically there would 
not be any difficulty for the passage of the particles. This 
reflection is partial and it becomes larger the closer the 
energy is closer to the height of the potential step. We 
can say that a particle with energy slightly positive is 
almost as confined as a particle inside the well. From 
this fact results the existence of almost bound states of 
positive energy known as quasi- stationary states or reso- 
nances. These resonances appear as peaks in the excita- 
tion function, a peak at a given energy meaning that the 
energy coincides with a given resonance of the nucleus. 

The existence of resonances can also be inferred from 
the properties of the wave function. We consider only 
elastic scattering, with the other channels closed. The 
external and internal wavefunctions are both sine func- 
tions, the first with wavenumber k = \/2mE/h and the 
second with K = ^2m(E — Vo)//i. If E is small and Vq 
is about —35 MeV, we have K » k. The internal and 
external parts should join at r = i? with continuous func- 
tion and derivatives. As the internal frequency is much 
larger than the external one, the internal amplitude is 
quite reduced. Only at the proximity of the situation in 
which the derivative is zero there is a perfect matching 
between both and the internal amplitude is identical to 
the external one. The energy for which this happens is 
exactly the energy of resonance. 

Resonances appear in the excitation function at rela- 
tively low energies, where the number of open channels 
is not very large and it is possible that to return to the 
entrance channel. To arrive at an expression of the cross 
section that describes a resonance, we rewrite [35] as. 



ui ~ exp{—iKr) -\- b exjp{iKr) , 



(r < R), (39) 



this time containing a second part which takes into con- 
sideration the part of the wave that returns. This second 
part allows the existence of resonant scattering, where 
the incident particle is re-emitted with the same energy 
that it entered, after forming the compound nucleus. The 
complex amplitude b is always smaller than one, because 
there are no creation of particles in the region < R in eq. 

m 

The second parenthesis in the denominator of eq. [26] 
is never zero, since the numerator forces // to be always 
negative. If for a certain energy fn vanishes, ar,o passes 
by a maximum in that energy. We can tentatively iden- 
tify these energies as being the energy of the resonances. 
Let us take the extreme case of a single resonance at the 
energy Ej^, that is, fn = ior E = Ej^. We can expand 
fji in a Taylor series in the neighborhood of a resonance, 
fniE) = {E-ER){dfR/dE)E=En + Keeping just the 



first term of the expansion and using 24 and 26 we get 



exp(2ifci?) - 1 



iVa 



Eu)+il 



(40) 



where we define: 



ttA^ 



r«(r-r„ 



[E - ErY 



(41) 



2kR 
{dfR/dE) 



and r = 



E=Ee 



2kR - 2fi 
{dfR/dE) 



E = Er, 



(42) 

The energy F that appears in [40] is the total width 
of the resonance, F = F^ -I- F^ -I- i.e., the sum of the 
widths for all the possible processes of decay of the nu- 
cleus, starting from the resonant state. F^ is the entrance 
channel width, and F — Fq, is the sum of the widths of all 
the exit channels except a. If we restrict the exit chan- 
nels to a single channel (3, or we denote /3 as the group 
of exit channels except a, eq. [41] is rewritten as 



7a, 13 



F„F 



iE-ERr + {^y 



(43) 



which is the usual way of presenting the Breit-Wigner 
formula, which describes the form of the cross section 
close to a resonance. Let us recall that eq. 43 refers to 
an incident particle of / = 0, without charge and without 
spin. If the spins of the incident and target particles 
are Sa and sa, respectively, and the incident beam is 
described by a single partial wave I ^ 0, one can show 
that the cross section of eq. [43] should be multiplied by 
the statistical factor g = (2/-|-l)/(2sa-|-l)(2s^-t-l) where 
/ is the quantum number of the total angular momentum 
I = Sa -I- Sa -I- 1 of the compound nucleus, g reduces, 
naturally, to the unit in the case of zero intrinsic and 
orbital angular momenta. 

If the exit channel is the same as the entrance chan- 
nel a, the cross section should be obtained from [40] and 



its dependence in energy is more complicated because in 
addition to the resonant scattering there is the poten- 



tial scattering, and the cross section 40 will contain an 
interference term between both. The presence of inter- 
ference results in a peculiar aspect of the scattering cross 
section, which differs from the simple form [43[ for the re- 
action cross section. This is seen in figure ]8|that shows 
the forms that a resonance can take in the scattering 
cross section. 

The region of energy where resonances show up can 
extend to 10 MeV in light nuclei but it ends well be- 
fore this in heavy nuclei. Starting from this limit the 
increase in level density with energy implies that the av- 
erage distance between levels is smaller than the width of 
the levels and individual resonances cannot be resolved 
experimentally. They form a continuum and this region 
is known as continuum region. The cross-section in this 
region fluctuates. 



D. Statistical theory of nuclear reactions 

The fluctuating behavior of low-energy nuclear reac- 
tions is due to the interference of the reaction amplitudes 
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FIG. 8; Differential cross section at four angles for the elastic 
scattering of protons off ^'^Al, in units of the Rutherford cross 
section, in the neighborhood of the 1579 keV resonance. For 
larger angles we have typical interference between the reso- 
nant scattering and the potential scattering [6]. 
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FIG. 9: Total cross section for neutrons incident on cad- 
mium, holmium and lead, showing an oscillatory behavior. 
The curves for cadmium and lead are results of calculations 
with the optical model [Tj. 



corresponding to the excitation of each of the overlap- 
ping states which vanish in the energy average of the 
cross-section since these amplitudes are complex func- 
tions with random modulus and phase. Calling crjv(c) 
the cross section for the formation of a compound nu- 
cleus in the entrance channel c, and using the reciprocity 
theorem which relates the cross-section acc' to the cross- 
section for the time-reversed process c' — >■ c, one gets 



acc' {Ec')dEc' = acN (c) 

{2Ic' + l)^ic'Ec' acN {c')uj{Uc')dUc' 
/(f (2/c + 1) l^cEc acN (c) ^{Uc)dUc 



,(44) 



where Ic is the angular momentum and /ic is the reduced 
mass in channel c. Ejectiles with energy in the range Ec' 
to Ec' + dEc' leave the residual nucleus with energy in 
the range Uc' to Uc' + dUc' where Uc' = Eqn — Be' — Ec' 
and Ec n and Be' are respectively the compound nucleus 
energy and the binding energy of the ejectile in the com- 
pound nucleus. Eq. |44] is the Weisskopf-Ewing formula 
for the angle- integrated cross-sections [TT] . To a good ap- 
proximation, the level density uj{U) cx exp(?7/T), so the 
ejectile spectrum given by the Weisskopf-Ewing theory is 
Maxwellian. It rises rapidly above the threshold energy, 
attains a maximum and then falls exponentially. 

The Weisskopf-Ewing theory simple to use, but it has 
the disadvantage that it does not explicitly consider the 
conservation of angular momentum and does not give the 
angular distribution of the emitted particles. This is pro- 
vided by the Hauser-Feshbach theory |12j . This theory 



takes into account the formation of the compound nu- 
cleus in states of different J and parity tt. Let us consider 
the case of a reaction leading from the initial channel c 
to a final channel c'. If there is no pre- equilibrium emis- 
sion, one may identify the compound nucleus formation 
cross-section (Jcn — I^Jtt'^cn with the optical model 
reaction cross-section = ^ X)/ (2Z + 1) T; , which, if 
the transmission coefficients T/ = 1 — KS*;)]^ do not de- 
pend on J. {Si) is the average value of the scattering 
amplitude over several overlapping resonances. 

Using the above assumptions, the reciprocity theorem, 
and following a similar derivation as with the Ewing- 
Weisskopf theory, the cross-section for transition c — > c' 
is given by the Hauser-Feshbach expression 



-y 



(2J + 1) Es.iTi{c)Es',i'Ti'{c') 



fc2 



(2i, + 1) (2/, + 1) EcEs.iTiic) 



(45) 



The compound nucleus states may be both of positive 
and negative parity. Since parity is conserved, in evalu- 
ating (45), one must take into account that the parity of 



compound nucleus states and the parity of the residual 
nucleus states may impose restrictions to the values of 
the emitted particle angular momentum. Thus, positive 
parity compound nucleus states decay to positive parity 
states of the residual only by even angular momenta and 
to negative parity residual nucleus states by odd angular 
momenta. 
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E. The optical model 



Expression 18 shows that if a particle is always ab- 
sorbed when it reaches the nucleus, the total cross section 
falls monotonically with the energy and grows linearly 
with . In figure [9] we see curves of total cross section 
of neutron scattering showing an oscillatory behavior for 
its energy dependence as well as for the mass dependence 
of the target. Their presence is mainly due to interfer- 
ence phenomena between the part of the incident beam 
that passes through the nucleus and the part that passes 
around it [8]. 

The basis of the optical model was established by Her- 
man Feshbach and collaborators in 1953 ^9^. In this 
model the interaction between the nuclei in a reaction 
is described by a potential U{r), with r being the dis- 
tance between the center of mass of the two nuclei. One 
replaces the complicated interaction that a nucleon has 
with the rest of the nucleus with a potential that acts on 
the nucleon. The potential U{r) includes a complex part 
that takes into account the absorption effects, i.e., the 
inelastic scattering. The nuclear scattering is treated in 
similar form as the scattering of light by a glass sphere 
and the name of the model derives of this analogy. 

In its most commonly used form, the optical potential 
is written as the sum: 

f/(r) = Unir) + Ui{r) + Uoir) + Us{r) + Uc{r), (46) 

which contains parameters that can vary with energy and 
masses of the nuclei and that should be chosen by an ad- 
justment to the experimental data. Obviously, the opti- 
cal potential U{r) will only make sense if these variations 
are small for neighboring masses or energies. 

The first part of |46) Unir) = -Vf{r,R,a), is real 
and represents a nuclear well with depth ~V, being 
multiplied by a Woods-Saxon form factor /(r, R, a) — 
[1 + exp{(r — i?)/a}]^ , where R is the radius of the nu- 
cleus and a measures the diffuseness of the potential, i.e., 
the width of the region where the function / is sensibly 
different from or 1. V, R and a are treated as adjustable 
parameters. 

The absorption effect or, in another way, the dis- 
appearance of particles from the elastic channel, is 
taken into account including the two following imag- 
inary parts Ui{r) — ~iWf{r,Ri,ai) and Uoir) — 
'iiaiWDdf{r,Rj,ai)/dr. An imaginary part produces 
absorption. It is easy to see this for the scattering 
problem of the square well: if an imaginary part is 
added to the well, U{r) — —Vq — iWq (r < tq) and 
t/(r) = (r > To), it appears in the value of 

K = [2m{E -h Vo + iWo)]^/V^- This wiU produce an 
exponentially decreasing internal wavefunction. Thus, it 
corresponds to an absorption of particles from the inci- 
dent beam. 

Ui is responsible for the absorption in the whole vol- 
ume of the nucleus, but C/^j, built from the derivative of 
the function /, acts specifically in the region close to the 
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FIG. 10: (Top) Angular distribution of the elastic scattering 
of 17 MeV protons on nuclei in the region Z = 26 — 30. (Bot- 
tom) The polarization value [48| in the scattering of 9.4 MeV 
protons on copper. The curves are, in both cases, obtained 
with fits from the optical model [14j . 



nuclear surface. These two parts have complementary 
goals: at low energies there are no available unoccupied 
states for nucleons inside the nucleus and the interactions 
are essentially at the surface. In this situation Uo{r) is 
important and Ui{r) can be ignored. On the other hand, 
at high energies the incident particle has larger penetra- 
tion and in this case the function Ui{r) is important. 
As with the shell model potential, used in nuclear struc- 
ture, a spin-orbit interaction term is added to the optical 
potential. This term, which is the fourth part of |46j is 
usually written in the form 

Us{r)^s-\Vs-^f[T,Rs,as). (47) 
r dr 

s is the spin operator and 1 the angular orbital momen- 
tum operator. As with Uoir), the part Us{r) is only 
important at the surface of the nucleus since it contains 
the derivative of the form factor /. The values oiVs, Rs 
and 05 must be adjusted by experiment. 
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The presence of the term Us is necessary to describe 
the effect of polarization. Through experiences of double 
scattering it can be verified that proton or neutron beams 
suffer strong polarization at certain angles. This means 
that the quantity 



P = 



Nc-Nb 



(48) 



where Nc is the number of nucleons in the beam with 
spin up and Nf, with spin down, has a value significantly 
different from zero at these angles. With the inclusion of 
Us, the optical model is able to reproduce in many cases 
the experimental values for the polarization |48[ 

Finally, a term corresponding to the Coulomb potential 
is added to |46] whenever the scattering involves charged 
particles. It has the form 



Ucir) = 



(r < R,) (49) 
(r>i?,), (50) 



where is assumed that the nucleus is a homogeneously 
charged sphere of radius equal to the Coulomb barrier 
radius Rc, which separates the regions of nuclear and 
Coulomb forces. 

Figure [TO] exhibits the result of the application of eq. 
[46|to the elastic scattering of 17 MeV protons on several 
light nuclei. The angular distribution is very well repro- 
duced by the model, which also reproduces correctly the 
polarization [46] for copper as a function of the scattering 
angle. 

The optical model has a limited set of adjustable pa- 
rameters and is not capable to describe abrupt varia- 
tions in the cross sections, as it happens for isolated res- 
onances. However, it can do a good description of the 
cross sections in the presence of the oscillations of large 
width in the continuous region, as it treats these as wave 
phenomena. 



IV. DIRECT REACTIONS 

Direct reactions becomes more probable as one in- 
creases the energy of the incident particle: the wave- 
length associated to the particle decreases and localized 
areas of the nucleus can be "probed" by the projectile. In 
this context, peripheral reactions, where only a few nucle- 
ons of the surface participate become important. Direct 
reactions happen during a time of the order of 10~^^s. 
Reactions with formation of compound nuclei can be up 
to six orders of magnitude slower. A reaction type at a 
given energy is not necessarily exclusive; the same final 
products can be obtained, part of the events in a direct 
way, other parts through the formation and decay of a 
compound nucleus. 

There are two characteristic types of direct reactions. 
In the first the incident particle scatters inelastically and 
the transferred energy is used to excite a collective mode 
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FIG. 11: Angular distribution of the reaction ^^P(d,n)^^S, 
with the transfer of a proton to several states of '^^S. The 
curves are results of DWBA calculations for the indicated I 
values [T5] . 



of the nucleus. Rotational and vibrational bands can be 
studied in this way. The second type involves a modifi- 
cation in the nuclear composition. Examples are transfer 
of nucleons, as pick-up and stripping reactions. An im- 
portant reaction of the latter kind is a knock-out reaction 
where the incident particle knocks out a particle of the 
target nucleus and continues in its path, resulting in three 
reaction products. Reactions with nucleon exchange can 
also be used to excite collective states. An example is 
a pick-up reaction where a projectile captures a neutron 
from a deformed target and the product nucleus is in an 
excited state belonging to a rotational band. 

Direct reactions exhibit a peculiar form of angular 
distribution, which allows us to extract information on 
the reaction mechanism with the employment of sim- 
ple models. Typical examples are the stripping reactions 
(d,n) and (d,p), where the angular distribution of the re- 
maining nucleon presents a forward prominent peak and 
smaller peaks at larger angles, with the characteristic as- 
pect of a diffraction figure. 

Figure |TT] shows experimental results for the reaction 
'^^Pfd.n)'^^ [T5]: angular distributions of the detected 
neutrons corresponding to each energy level of '^^S are 
exhibited for several angular momentum. We see that 
the behavior of the cross sections is in agreement with 
qualitative predictions: the curves exhibit a first peak at 
a value of 9 that grows with I. Other smaller peaks occur 
as 9 increases. The increase of 9i with I is an important 
characteristic that can be used to identify the value of the 
transferred momentum in a given angular distribution. 

We consider an initial quantum state with particles 
of mass nia hitting a target A of mass tua- The final 
quantum state are particles of mass mt, moving away from 
nucleus B of mass ms- The differential cross section for 
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the process is given by 

da _ mgmbkb , ,2 

involving the matrix element 



(51) 



(52) 



^a,^b, ^A, are the internal wavefunctions of the nu- 
clei a, b, A and B. 5'q,^/3 are the wavefunctions of the 
relative motion in the entrance channel a and in the exit 
channel /3. The integration volume dr spans the coor- 
dinates of all particles. V is the perturbation potential 
that causes the "transition" from the entrance to the exit 
channel. 

The use of plane waves for 'i'a and 'i'p in eq. 52 is 
known as first Born approximation. With it we can ar- 
rive to an approximate expression for the behavior of the 
differential cross section. As the nuclear forces are of 
short range, one can restrict the integral [52] to regions 
where = = r. This leads to 



Vf,^ y"drexp(ik-r)|y"^'^*^^*,«'AdT' 



(53) 



where k = ka— k^. The global variables in dr were sep- 
arated into variables dr and dr' . 

Expanding the plane wave in a Legendre polynomial 
series, we obtain 



00 „ 
l^/, ?^^^i'(2/ + l) / j,{kr)Pi{cos0)F{r)dr, 
1=0 •' 



(54) 



where F(v) = J **iJ/^T/vl/^vl/^dr' contains all the inter- 
nal properties and is known as as the form factor of the 
reaction. The action of V is restricted to the surface of 
the nucleus: outside the nucleus the action of V is lim- 
ited by the short range of the nuclear forces and inside 
the nucleus there is a strong deviation to the absorption 
channel. Expression 54 becomes Vfi = J2'iZo'^iji{kR)y 



becomes Vji 

where the coefficients q contain information on the form 
factor F{r). The index I can be identified as the angular 
momentum transferred and, for a reaction that involves 
a single value of I, we can write for the differential cross 
section; 



da 

dn 



(X \ji{kR)\' 



where the dependence in 9 is contained in 



kp — 2kakp cos 6. 



(55) 



(10.82) 



We have an oscillatory behavior for the angular distribu- 
tion, the maxima separated by tt from each other in the 
axis kR. 

The Born approximation with plane waves predicts for 
certain cases the correct place of the first peaks in the 



angular distribution but without reproducing correctly 
the intensities. A considerable progress can be done in 
the perturbative calculations if, instead of plane waves 
in |52[ we use distorted waves that contain, besides the 
plane wave, the part dispersed elastically by the optical 
potential of the target. The Born approximation with 
distorted waves, or DWBA {distorted wave Born approx- 
imation) , became a largely employed tool in the analysis 
of experimental results of direct reactions. With it one 
can try to extract with a certain reliability the value of 
the angular momentum I transferred to the nucleus in 
a stripping or pick-up reaction. An example of this is 
the already mentioned stripping reaction, '^^P(d,n)^^S, 
for deuterons of 25 MeV. For the energy levels shown in 
figure [11] the assignments of the value of I for the level is, 
in most cases, univocal. 

The angular momentum I transferred in a direct re- 
action generally modifies the value of the total angular 
momentum of the nucleus. If Ji is the spin of the target 
nucleus, the spin Jf of the product nucleus is limited to 
the values 



< Jf < J, + l + 



1 



(56) 



and the initial an d fi nal parities obey the relationship 

TT.nj = i-iy. 



Eq. 



56 



allows, with the knowledge of the 
target nucleus and of the transferred angular momen- 
tum, the determination of the parity of the product state 
formed and is a tool for the determination of its spin. 

The knowledge of the transferred angular momentum 
value in a direct reaction opens the possibility to test 
the predictions of the shell model for the structure of nu- 
clei. In a direct reaction one assumes that the nucleon 
is located in an orbit of the nucleus with the same an- 
gular momentum as the transferred momentum in the 
reaction. In almost all cases studied with direct reac- 
tions the value of the assigned / corresponds exactly to 
the predicted by the shell model. We know, however, 
that the real situation is more complicated, due to the 
presence of the residual interactions that give place to 
configuration mixing. As result, the cross section for the 
formation of a state i of the product nucleus is related 
to that calculated with DWBA for the formation from a 
single-particle state by 



(-) 
\ dfl J exp 



2Ji3 + l fda 



2Ja + 1 ^KdQJBWBA' 



(57) 



where the spectroscopic factor Sij measures the weight 
of the configuration j used in the DWBA calculation, in 
the final state «, with the sum-rule 



(58) 



The sum [58] embraces all the nucleons in the configura- 
tions i of the product nucleus. The statistical weight 
{2Jb + 1)/(2Ja + 1) that appears in the DWBA calcu- 
lation involving the angular momentum of the target nu- 



cleus J A and final nucleus Jb, is explicitly given in 57 
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V. HEAVY ION REACTIONS 

A. Types and properties 

Heavy ion reactions (with A > 4) can be separated 
into 3 major categories. 

1) Due to their large charge, two heavy nuclei feel a 
strong mutual Coulomb repulsion. To produce a nuclear 
reaction the projectile needs enough energy to overcome 
the Coulomb barrier. For a very heavy target, as ^'^^U, it 
is necessary about 5 MeV per nucleon. Then the wave- 
length of the projectile is small compared with the di- 
mensions of the nuclei and classical and semi-classical 
methods become useful in the description of the reaction. 

2) The projectile carries a large amount of angular mo- 
mentum and a good part of it can be transferred to the 
target in the reaction. Rotational bands with several 
dozens of units of angular momentum can be created. In 
fact, heavy ion reactions are the best suited to feed high 
spin levels. 

3) Direct reactions and formation of compound nucleus 
are also common processes in reactions with heavy ions. 
But some peculiarities of these are not found in reactions 
with projectile nucleons. One of these processes can be 
understood as intermediate between a direct reaction and 
the formation of a compound nucleus. Fusion does not 
occur but projectile and target pass a relatively long time 
under the mutual action of the nuclear forces. Nuclear 
matter is exchanged between both and there is a strong 
heating of the two nuclei, with a large transfer of kinetic 
energy to the internal degrees of freedom. These are the 
deep inelastic collisions. 

The kind of process that prevails depends upon the 
distance of closest approach d between the projectile 
and target. If this distance is sufficiently large only 
the long range Coulomb interaction acts and, for a clas- 
sical hyperbolic trajectory, d is related to the impact 
parameter b and to the energy E of the projectile by 

d = a/2 + (a/2)^ + where a is the distance of 

closest approach in a head-on collision. It is this is re- 
lated to Ehy a = ZiZ2e^ /AneoE. 

Experimentally, the variable under control is the en- 
ergy E of the projectile and, for E sufficiently large, d 
can be small enough to enter the range of nuclear forces. 
Collisions near this limit are called grazing collisions and 
are characterized by values of bgraz and dgraz ■ Assuming 
that there is always reaction when b < bgraz y the reac- 
tion cross section ar can be determined geometrically by 
CTr = '^b'^raz- The experimental determination of ar al- 
lows to establish the value dras = 0.5 + 1.36(ylJ^'^ 
showing that the distance of grazing collision is some- 
what larger than that deduced from two touching spheres 
(1.36 fm > ro = 1.2 fm). 

When the impact parameter is close to bgraz one ex- 
pects nuclear reactions of short duration, without the 
contribution of the compound nucleus formation. Such 
reactions are elastic and inelastic scattering and transfer 



of few nucleons. When the incident energy is sufficiently 
high, small values of b can lead to the projectile pene- 
trating the target. Depending on the energy and on the 
involved masses, the reaction can end in one of the pro- 
cesses below: 

a ) Fusion - is the preferred process when one has light 
nuclei and low energy. There is the formation of a highly 
excited compound nucleus that decays by evaporation 
of particles and 7-radiation emission, leading to a cold 
residual nucleus. If the energy in the CM is close to the 
Coulomb barrier energy the cross section of compound 
nucleus formation starting from two nuclei is practically 
equal to the reaction cross section. 

b ) Fission - When the compound nucleus is heavy the 
fission process competes strongly with the evaporation of 
particles in each stage of the evaporation process. A very 
heavy compound nucleus with a large excitation energy 
has a very small probability of arriving to a cold residual 
nucleus without fission at some stage of the de-excitation. 
The role of the angular momentum I transmitted to the 
target nucleus is also essential. The fission barrier de- 
creases with the increase of I and for a critical value Iq^^^^ 
the barrier ceases to exist. A nucleus with angular mo- 
mentum greater than suffers immediate fission and 
this is also a limiting factor in the production of super- 
heavy elements. 

c) Deep inelastic collision (DIG) - is a phenomenon 
characteristic of reactions involving very heavy nuclei 
[A > 40) and with an incident energy of 1 MeV to 
3 MeV above the Coulomb barrier. In DIC the projectile 
and the target spend some time under mutual action, ex- 
changing masses and energy but without arriving to the 
formation of a compound nucleus. The projectile escapes 
after transferring part of its energy and angular momen- 
tum to its internal degrees of freedom and to the target, 
with values reaching 100 MeV and 50h. 

One of the most interesting aspects of DIC is the cor- 
relation between the energy dissipated in the collision 
and the scattering angle in the center of mass. Let us 
look at figure [T2|top). The trajectory 1 shows the pro- 
jectile with an impact parameter that leads it out of the 
range of nuclear forces. The Coulomb scattering angle 
will become larger as the impact parameter decreases. 
In a graph as in figure 12 bottom), where one plots the 



final energy against the scattering angle, trajectories of 
type 1 are located in the upper branch, where there is 
no dissipation and the initial kinetic energy stays unaf- 
fected. This upper branch has a maximum value for the 
scattering angle. At a smaller impact parameter the nu- 
clear attractive force begins to act and, with that the 
effects of dissipation of DIC. A given Coulomb scattering 
angle 9 can also be reached by the combination of the 
nuclear and Coulomb forces. Only that now there is loss 
of energy and the events are located in the branch 2 of 
figure [l2|b) . There still is no one-to-one correspondence 
between angle and energy because an infinite number of 
trajectories can lead to the same angle 9. The branch 
2 should be understood as a line of maxima in a three- 
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FIG. 12: (Top) Reactions with different impact parameters 
leading to the same scattering angle Q. (Bottom) Lines of 
maxima in topographical diagrams of the final energy against 
the scattering angle. 



dimensional representation (called Wilczynski diagram) 
where the axis perpendicular to the paper is proportional 
to the cross section cPa/dEdO. The same angle can also 
be obtained by the trajectory 3, with a longer interaction 
time between the nuclei and a larger dissipation. As now 
the projectile is deflected towards the nucleus, the scat- 
tering angle would be formally —9 but, as there is no 
experimental distinction between and —0, these events 
appear as an independent ramification in the lower part 
of the diagram. 



B. Superheavy elements 

The heaviest element found in the nature is ^gfU. It 
is radioactive, but it survived since its formation in su- 
pernovae explosions because it has a decay half-life of 
the order of the age of the Earth. Elements with larger 
atomic number (transuranic) have shorter half-lives and 
have disappeared. They are created artificially through 
nuclear reactions using heavy elements as target. Ini- 
tially, the projectiles used were light particles: protons, 
deuterons, a-particles and neutrons. The use of neutron 
is justified because the /3^-emission of the compound nu- 
cleus increases the value of Z and it was in this way jl9] 
that the first transuranic element, the neptunium, was 
obtained: tl+ 238u ^ ^sgy ^ "^^^^^ + p- . 

Reactions with light particles can produce isotopes up 
to mendelevium {Z = 101), but it is not possible to go 
beyond that; the half- lives for a-emission or spontaneous 



fission become extremely short, turning impracticable the 
preparation of a target. The alternative is to place a 
heavy element under the flux of very intense neutrons. 
This can be done using special reactors or using the rest 
material of nuclear explosions. The elements einsteinium 
{Z = 99) and fermium [Z = 100) were discovered in this 
way in 1955 but the increasing competition that the beta 
decay has with alpha decay and with spontaneous fission 
prevents this method to be used for larger Z. 

Starting in 1955 heavy ion accelerators began to deliver 
beams with high enough intensity and energy to compete 
in the production of transuranic isotopes. The first posi- 
tive result was the production of two californium isotopes 
{Z — 98) in the fusion of carbon and uranium nuclei: 

+ 2|«u ^ ^gtcf+en, 

'eC ^ '9'|Cf+4n. (59) 

This opened the possibility of reaching directly the nu- 
cleus one wants to create from the fusion of two smaller 
nuclei. The difficulty of such task is that the cross sec- 
tions for the production of heavy isotopes are extremely 
low. As example, the reaction |^Ti -I- ^°|Pb lolRf + 
n, which produces the element rutherfordium, has a cross 
section of only 5 nb. A small increase in the charges re- 
duces drastically this value: the cross section for the fu- 
sion reaction ^|Fe -I- ^gfPb -J> lolHs -I- n is 4 pb. As 
comparison, the typical cross sections of DIG for heavy 
nuclei are in the range 1-2 b. 

In spite of the experimental refinement that these low 
cross sections demand, one is able to produce isotopes 
with charge as heavy a.s Z — 118. The understand- 
ing of the mechanisms that lead to fusion is, however, 
not fully understood. According to the traditionally ac- 
cepted model, the fusion of two nuclei proceeds in two 
stages: the formation of a compound nucleus and the 
de-excitation of the compound nucleus by evaporation of 
particles, preferentially neutrons. The difficulties for the 
materialization of the process in very heavy nuclei reside 
in both stages pO] . 

VI. ELECTROMAGNETIC PROBES 
A. Coulomb Excitation 

Coulomb excitation is a inelastic scattering process in 
which a nucleus excites another nucleus with its elec- 
tromagnetic field V. This field can be decomposed in 
terms of a series of multipoles, e.g. El, E2, Ml, 
which carry well defined angular momenta and parities. 
At low bombarding energies E2 (quadrupole) excitations 
are more common, while at higher energies El excita- 
tions dominate. As an example of a low energy reaction, 
let us consider the excitation of a quadrupole state in 
a head-on collision below the Coulomb barrier, i.e. for 
a situation in which the projectile deccelerates as it ap- 
proximates the target and stops before reaching the range 
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of the nuclear interaction, reaccelerating backwards after 
that. The differential cross section is given by the prod- 
uct of the Rutherford differential cross section at 180° 
and the probability of transition of the target from state 
i to state /, along the trajectory, measured by the square 
of 



da 



?=180° 



d(7R 



e=i80° 



X |a,/|^ 



(60) 



The square of / measures the transition probability 
from i to f and this probability should be integrated along 
the trajectory. A simple calculation for a^/ can be done 
in the case of the excitation from the ground state J = 
of a deformed nucleus to an excited state with J = 2. The 
perturbation V comes, in this case, from the interaction 
of the projectile charge Zpe with the quadrupole moment 
of the target nucleus. This quadrupole moment works as 
an operator that acts between the initial and final states. 



V = 



(61) 



where r is the projectile-target separation distance and 



*}.(3z2 



(62) 



where the sum extends to all protons at the positions 
ri = (xijiji^Zi) in the target. The amplitude is then 
given by 
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(63) 



Using eq. [BOfor the Rutherford differential cross sec- 
tion at 6* = 180°, we obtain 
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(64) 



which is an expression that is independent of the charge 
of the projectile. It is, on the other hand, proportional 
to the mass of the projectile, indicating that heavy ions 
are more effective for Coulomb excitation. 

The quadrupole moment operator Qif uses the wave- 
functions and of the initial and final states. If 
those two wavefunctions are similar, as is the case of an 
excitation to the first level of a rotational band, the op- 
erator Qif can be replaced by the intrinsic quadrupole 
moment Q. The expression translates, in this way, the 
possibility to evaluate the quadrupole moment from a 
measurement of the cross section. This has been indeed 
a major spectroscopic tool for determining quadrupole 
momenta (and transitions) along the nuclear chart. 



B. Photonuclear reactions and giant resonances 

A photonuclear reaction is a reaction resulting from 
the interaction of the electromagnetic radiation with a 
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FIG. 13: Giant resonance in the absorption of photons by 
^20Sn in]. 



nucleus. Therefore, one can access information which are 
complementary to Coulomb excitation (and vice- versa). 

When the energy of the photon is located above the 
separation energy of a nucleon, the cross section of photo- 
absorption reveals the presence of characteristic sharp 
resonances. But when the incident energy reaches the 
range of 15-25 MeV, a new behavior appears in the cross 
section, with the presence of a wide and large peak, called 
giant electric dipole resonance {El excitation). Figure 
[13] exhibits the excitation function of photoabsorption of 
^^Sn at photon energies around the electric dipole giant 
resonance at 15 MeV. 

The giant resonance occurs in nuclei along the whole 
periodic table, with the resonance energy decreasing was 
Egdr - 80/^1/3 for A > 20. Their widths are almost all 
in the range between 3.5 MeV and 5 MeV. In few cases 
they can reach 7 MeV. They are a collective excitation, 
i.e., an excitation involving many nucleons at once. The 
time-varying electric field of the photon is very effective 
in inducing collective oscillations of protons against neu- 
trons. 

The giant electric dipole (GDR) resonance arises from 
an excitation that transfers by one unit of angular mo- 
mentum to the nucleus {Al = 1). If the nucleus is even- 
even it is taken to a 1^ state. The transition also changes 
the isospin of 1 unit (AT = 1) and, due to that, it is 
also called isovector resonance. The photon can excite 
less effectively, but still with appreciable cross sections, 
giant isoscalar resonances, with AT = 0. But electric 
quadrupole (AZ — 2) and electric monopole [Al = 0) 
are observed mostly in reactions with charged particles. 
In a giant electric quadrupole resonance the nucleus os- 
cillates between an spherical (assuming that this is the 
form of the ground state) to an ellipsoidal form. If pro- 
tons and neutrons vibrate in phase, we have an isoscalar 
resonance (AT = 0) and if they are oscillate in opposite 
phase we have an isovector resonance (AT — 1). The 
giant monopole resonance is a very special case of nu- 
clear excitation where the nucleus contracts and expands 
radially, maintaining its original form but changing its 
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volume. It is also called breathing mode. It can also hap- 
pen in isoscalar and isovector forms. It is an important 
way to study the compressibility of nuclear matter. 

Besides electric giant resonances, associated to a varia- 
tion in the form of the nucleus, magnetic giant resonances 
exist, involving spin vibrations. In these, nucleons with 
spin up move out of phase with nucleons with spin down. 
The number of nucleons involved in the process cannot 
be very large because it is limited by the Pauli principle. 
The magnetic resonances can also separate into isoscalar 
resonances, where protons and neutrons of same spin 
vibrate against protons and neutrons of opposite spin, 
and isovector, where protons with spin up and neutrons 
with spin down vibrate against their corresponding part- 
ners with opposite spins. The last cases, are effectively 
probed in charge-exchange reactions, i.e. when the pro- 
jectile charge changes to Z-1 or Z+1, as in (d,p) and (d,n). 
They are known as giant Gamow-Teller resonances. 

Giant resonances can also be produced in an excited 
nucleus, a case known as the Brink- Axel hypothesis 
A special case is when two giant resonances are excited in 
the same nucleus. The double giant dipole resonance (or 
multiphonon giant resonance was observed for the first 
time in reactions with double charge exchange induced 
by pious in the '^^S [2B1. But Coulomb excitation proved 
to be their best probe, as shown theoretically in refs. 
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where q"^ — {Pf — Pi) ~ {Ef — Ei) is the 4-momentum 
transfer squared, where (p/) and Ei [Ef) are the elec- 
tron incoming (outgoing) momentum and energy. Ao(g^) 
and Bo{q^), are functions of q^, are the form factors as- 
sociated with the charge and magnetization distribution 
of the nucleus respectively. For a spin-0 nucleus: 



Aoiq^ 



1 + /4Af2 



and Boiq^)^0, (66) 



where Gsiq^) is the Sachs charge form factor, and the 
factor 1 -I- q^/4M^ is the kinematical recoil correction. 
For g 0, ss q2, where q is the three-momentum 
transfer, and to leading order in powers of q: 



GE{q') 



Ge{<1^) = J PT{v)e'^-'d\. 



(67) 



C. Electron scattering 

Electron scattering experiments have provided a rich 
database on most stable or long-lived nuclei. Electrons 
are structureless point-like objects that only interact elec- 
tromagnetically. Therefore, electron scattering avoids 
the complexity of the strong interaction between the pro- 
jectile and the target, and provides clean information 
about the charge distribution in the nucleus. 

Considering relativity and the spin of the electrons, the 
differential cross section for elastic scattering by spinless 
point-like nuclei can be expressed by the Mott scattering 
formula (for simplicity, we use h = c = 1 units): 
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Apl sin-* {9/2) [I + {2po/M) sin^ {9/2)] 

(65) 

where Zq^e and M are the charge and mass of the tar- 
get nucleus respectively, and po is the momentum of the 
incoming electron. Since the nucleus is not a point-like 
particle, the formula is modified by adding the nuclear 
electric and magnetic form factors, which contain infor- 
mation of the charge and magnetization distributions in- 
side the nucleus. 

The Rosenbluth formula accounts for the electron spin 
and explicitly expresses the cross section for arbitrary 



One thus sees that electron scattering is related to the 
Fourier transform of the nuclear charge density, Pt{^)- 
If the function Ge{<^) is mapped for a sufficient num- 
ber of momentum transfers q^ the Fourier transform can 
be inverted and pt{y) can be mapped with precision. 
For heavy elements {A > 10), this technique shows that 
the charge radius is closely proportional to 1.12 A^^^ fm, 
where A is the nuclear mass number. This indicates the 
density saturation for the nuclear matter. 

When q2 = 0, inelastic electron scattering probes 
the same multipolarity transitions {El, E2, ■ ■ ■) as in 
Coulomb excitation or with real photons. However, 
in contrast to the later probes, electron scattering also 
provides information on nuclear excitations for cases in 
which q2 7^ 0. This is useful in many aspects. For exam- 
ple, at high momentum transfers, the momentum transfer 
q can be absorbed by a nucleon in the nucleus which ac- 
quires an energy of the order of q2/2m^, thus probing 
the effective mass of a nucleon in the nucleus. This is 
of relevance to understand the effects of nucleon-nucleon 
interactions in the nuclear environment. The g-region 
where this occurs is termed by quasi-free scattering. As 
electrons can penetrate the nuclei, they can also more ef- 
fectively probe monopole (or breathing mode) excitations 
of the nuclei, a case in which the nuclear mass distribu- 
tion vibrates radially. This is of importance for determin- 
ing the compressibility modulus of the nuclear matter. Fi- 
nally, in high-energy electron scattering the electron can 
penetrate deeply inside the nucleons and probe the spin 
and charge distributions of quarks and gluons inside the 
nucleon 1261. 
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VII. RELATIVISTIC NUCLEAR COLLISIONS 
A. Transport theories and equation of state 

As the bombarding energy in nucleus-nucleus collisions 
increases the structure aspects of the nuclei become less 
relevant. Except for the bulk properties of the nuclei 
(size and number of nucleons), the physics involved is 
primarily due to the individual, and sometimes collec- 
tive, hadronic collisions. Several theoretical tools are 
used to describe these reactions, and in particular we 
quote (a) time-dependent Hartree-Fock {TDHF), (b) anti- 
symmetrized molecular dynamics (AMD), (c) Boltzmann- 
Uehling-Uhlenbeck (BUU, etc. 

At intermediate energies of Eiab ^ 100 — 1000 
MeV/nucleon the nucleons and the products of their colli- 
sions can be described individually and their propagation 
can be described by semiclassical equations. One of such 
equations, and perhaps the most popular in such stud- 
ies, is the so-called Boltzmann-Uehling-Uhlenbeck (BUU) 
equation: 
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dl 
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Vr/ - VrC/ • Vr/ = 

J d^p2 J dn aNN {^) I vi - V2 1 

X {/{/2 [1 - /i] [1 - /2] - /1/2 [1 - ./(] [1 - m , (68) 

If dN is the number of particles in the volume ele- 
ment d^r and whose momenta fall in the momentum 
element d'^p at time t, then the distribution function 
/ (r, p, t) is given by dN — f (r, p, t) d^rd^p. Thus the 
BUU equation is an equation for the distribution func- 
tion / (r,p,t). To account for the effect of each particle 
interacting with all others, one introduces the concept of 
mean-field, U (r, p, t) . This mean-field exerts a force on 
each particle, given by —V^U (r, p, t) . Also, the momen- 
tum dependence of the potential introduces a dependence 
through the derivative —VpC/ (r, p, t) . 

Due to the nucleon-nucleon collisions, the distribution 
function within d'^rd'^p can also be modified by nucleons 
leaving (or entering) this volume. This is taken care by 
the collision term, i.e. the right-hand-side of the BUU 
equation, umn is the nucleon-nucleon differential scat- 
tering cross section, Vi and V2 are the velocities of two 
colliding nucleons. The first factor inside braces are for 
collisions populating the volume element and the second 
term for those depleting it. The factors (1 — /) account 
for Pauli blocking of final occupied states. The integrals 
average over scattering angle and over all collisions within 
d^rd^p. The BUU equation falls in the category of what 
one calls quantum transport theories. Hadronic trans- 
port theories have been quite successful in applications, 
describing a multitude of measured particle spectra. 

Eq. (68) needs as basic ingredients the mean field U 



and the cross section aNN- Because these two quantities 
are related to each other, one should in principle derive 
them in a self-consistent microscopic approach, as in the 
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FIG. 14: - Contour plots of baryon density in the reaction 
plane in Au -I- Au collisions at 400 MeV/nucleon. The dis- 
played contour lines are for the densities pj pa — 0.1, 0.5, 1, 
1.5, and 2 m\. 



Brueckner theory. However, in practice the simulations 
are often done with a phenomenological mean field and 
free nuclear cross sections. 

An important ingredient in the transport theory calcu- 
lations is the compressibility K of nuclear matter, which 
refers to the second derivative of the compressional en- 
ergy E with respect to the density: 



K = 9p' 



5p2 



(69) 



This is an important quantity, e.g., for nuclear astro- 
physics. Supernova models might or not lead to explo- 
sions depending on the value of K. The central collisions 
of heavy nuclei are one of the few probes of this quantity 
in the laboratory. The dependence of the calculations on 
K follow from the dependence of the mean field poten- 
tial U [U ^ E/A+ kinetic energy terms) on the particle 
density p. A typical parametrization for U is the Skyrme 
parametrizationt/ = apj p^ -['b{p/ p^Y ■ 

The output of eq. (68) is the distribution function 
/(r,p,t), which allows one to calculate many proper- 
ties of heavy-ion collisions. Let us quote collective flows, 
proton and neutron production rates, (sub-threshold and 
above threshold) pion and kaon yields, etc. Combining 
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eq. (681 with a phase-space coalescence models one can 



also calculate such quantities as exclusive flows and in- 
termediate fragment formations. 

The dynamics of the central high-energy reactions can 
be broken down into several stages. Baryon-density 
contour-plots are shown in fig. 14 for 400 MeV/nucleon 
Au-|-Au collisions at & = 0, which will serve to illustrate 
our points. 

Following an initial interpenetration of projectile and 
target densities, the NN collisions begin to thermalize 
matter in the overlap region making the momentum dis- 
tribution there centered at zero momentum in the CM. 
The density in the overlap region rises above normal and 
a disk of excited and compressed matter forms at the 
center of the system. More and more matter dives into 
the region with compressed matter that begins to expand 
in transverse directions. At late stages, when the whole 
matter is excited, transverse expansion predominates. 
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B. Kinematics 

In relativistic nucleus- nucleus collisions other defini- 
tions are best suited for discussing relations between en- 
ergies, momenta and angles. The rapidity is a variable 
frequently used to describe the behavior of particles in 
inclusively measured reactions. It is defined by 



y 



2 \E-p« 



(70) 



which corresponds to tanhy = -g-, where y is the rapid- 
ity, p|| is the longitudinal momentum along the direction 
of the incident particle, E is the energy, both defined for 
a given particle. The accessible range of rapidities for a 
given reaction is determined by the available center-of- 
mass energy and all participating particles' rest masses. 
One usually gives the limit for the incident particle, elas- 
tically scattered at zero angle: 
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(71) 



where /3 = v is the velocity and all variables referring to 
the through-going particle given in the desired frame of 
reference (e.g. in the center of mass). 

Note that dy/dp\\ = 1/E. A Lorentz boost (3 along 
the direction of the incident particle adds a constant. 
In (7 -1-7/3), to the rapidity. Rapidity differences, there- 
fore, are invariant to a Lorentz boost. Statistical particle 
distributions are flat in y for many physics production 
models. Frequently, the simpler variable pseudorapidity 
1] is used instead of rapidity (and sloppy language mixes 
up the two variables). 

The pseudorapidity is a handy variable to approximate 
the rapidity if the mass and momentum of a particle are 
not known. It is an angular variable defined by 
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(72) 



FIG. 15: - As water comes in different phases (solid, liquid, 
gas) , so nuclear matter can come in its normal hadronic form 
and at sufficiently high temperature and density, in the form 
of a deconfined state of quarks and gluons. The diagram 
shows how nuclear matter should behave as a function of den- 
tistry and temperature. 



whose inverse function isO — 2 arctan (e^''), where 9 is 
the angle between the particle being considered and the 
undeflected beam. 77 is the same as the rapidity y if 
one sets f3 — 1 (or ni — 0). Statistical distributions 
plotted in pseudorapidity rather than rapidity undergo 
transformations that have to be estimated by using a 
kinematic model for the interaction. 



C. The quark-gluon plasma (QGP) 

The primary motivation for studying ultra-relativistic 
heavy ion collisions is to gain an understanding of the 
equation of state of nuclear, hadronic and partonic mat- 
ter, commonly referred to as nuclear matter. Displayed 
in fig. [15] is a schematic phase diagram of nuclear mat- 
ter. The behavior of nuclear matter as a function of 
temperature and density (or pressure), shown in fig. 15 
is governed by its equation of state. 

Conventional nuclear physics is concerned primarily 
with the lower left portion of the diagram at low temper- 
atures and near normal nuclear matter density. Here nor- 
mal nuclei exist and at low excitation a liquid-gas phase 
transition is expected to occur. This is the focus of ex- 
perimental studies using low energy heavy ions. At some- 
what higher excitation, nucleons are excited into baryonic 
resonance states, along with accompanying particle pro- 
duction and hadronic resonance formation. In relativistic 
heavy ion collisions, such excitation is expected to create 
hadronic resonance matter. 

We now briefly discuss the QGP signatures in nucleus- 
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FIG. 16: - Formation and evolution ol J/^ particles in rel- 
ativistic heavy ion collisions. After the formation, the J/3's 
are dissociated in the plasma due to color screening. The end 
effect is a smaller number J/^s of than expected from pure 
hadron-hadron multiple collisions. 



nucleus collisions. For more details see, e.g., ref. [251l29j . 

One group of such signatures can be classified as ther- 
modynamic variables. This class involves determination 
of the energy density e, pressure P , and entropy density 
s of the interacting system as a function of the temper- 
ature T and the baryochemical potential ■ Experi- 
mental observables can be identified with these variables 
and thus their relative behavior can be determined. If a 
phase transition to QGP occurs, a rapid rise in the effec- 
tive number of degrees of freedom, expressed by e/T^ or 
s/T^, should be observed over a small range of T. The 
variables T, s, and e, can be identified with the average 
transverse momentum (pt), the hadron rapidity density 
dN/dy, and the transverse energy density dEx/dy, re- 
spectively. The transverse energy produced in the inter- 
action is Et = ^j-EiSin^j , where E^ and 9i are the 
kinetic energies of the ejectiles and the emission angles. 

Electromagnetic (EM) probes, such as photons and 
leptons, provide information on the various stages of the 
interaction without modification by final state interac- 
tions. These probes may provide a measure of the ther- 
mal radiation from a QGP, if a region of photon energy, 
or equivalently lepton pair invariant mass, can be iso- 
lated for emission from a QGP relative to other processes. 
However, the yields for EM probes are small relative to 
background processes, which are primarily EM decays of 
hadrons and resonances. Lepton pairs from the QGP 
are expected to be identifiable in the 1-10 GeV invariant 
mass range. The widths and positions of the p, oj, and 
(j) peaks in the lepton pair invariant mass spectrum are 
expected to be sensitive to medium-induced changes of 
the hadronic mass spectrum 

The production of J/^* particles in a quark-gluon 
plasma is predicted to be suppressed (see fig. 



16). This 



is a result of the Debye screening of a cc pair, initially 
formed in the QGP by fusion of two incident gluons. Less 



tightly bound excited states of the cc system, such as ^' 
and Xc J E^re more easily dissociated and will be sup- 
pressed even more than the J/ 4". 

A long-standing prediction for a signature of QGP for- 
mation is the enhancement of strange hadrons. The 
production of strange hadrons relative to nonstrange 
hadrons is suppressed in hadronic reactions. This sup- 
pression increases with increasing strangeness content of 
the hadron. In a QGP the strange quark content is 
rapidly saturated by ss pair production in gluon-gluon 
reactions, resulting in an enhancement in the produc- 
tion of strange hadrons. Thus, multi-strange baryons 
and strange antibaryons are predicted to be strongly en- 
hanced when a QGP is formed. 

The connection between energy loss of a quark and 
the color-dielectric polarizability of the medium can be 
established in analogy with the theory of electromag- 
netic energy loss. Although radiation is a very effi- 
cient energy loss mechanism for relativistic particles, it is 
strongly suppressed in a dense medium by the Landau- 
Pomeranchuk effect |28J. Adding the two contributions, 
the stopping power of a quark-gluon plasma is predicted 
to be higher than that of hadronic matter. A quark 
or gluon jet propagating through a dense medium will 
not only loose energy but will also be deflected. This 
effect destroys the coplanarity of the two jets from a 
hard parton-parton scattering with the incident beam 
axis. The angular deflection of the jets also results in 
an azimuthal asymmetry. The presence of a quark-gluon 
plasma is also predicted to attenuate the emission of jet 
pairs opposite to the trigger jet. 

All of the above quark-gluon plasma signatures have 
been studied to maturity at the CERN (European Or- 
ganization for Nuclear Research) - SPS (Super Proton 
Synchrotron) and the BNL (Brookhaven National Lab- 
oratory) - RHIC (Relativistic Heavy Ion Collider) facil- 
ities. For more details, the following references discuss: 
(a) direct photons [301111], (b) dileptons [12, (c) J/^- 
production !33j and (d) quark in-medium attenuation, 
also known as high-py quenching or jet attenuation [34] . 



VIII. NUCLEAR REACTIONS IN STARS 
A. Hydrogen and CNO cycles 

Energy production in stars is a well known process. 
The initial energy which ignites the process arises from 
the gravitational contraction of a mass of gas. The con- 
traction increases the pressure, temperature, and density, 
at the center of the star until values able to start the 
thermonuclear reactions, initiating the star lifetime. The 
energy liberated in these reactions yield a pressure in the 
plasma, which opposes compression due to gravitation. 
Thus, an equilibrium is reached for the energy which is 
produced, the energy which is liberated by radiation, the 
temperature, and the pressure. 

The Sun is a star in its initial phase of evolution. The 
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FIG. 17: The p-p chain reaction (p-p cycle). The percentage 
for the several branches are calculated in the center of the 
Sun ED. 
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FIG. 18: The CNO cycle, (courtesy of Frank Timmes). 
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temperature in its surface is 6000° C, while in its interior 
the temperature reaches 1.5x 10'' K, with a pressure given 
by 6x 10^^ atm and density 150 g/cm'^. The present mass 
of the Sun is Mq = 2 x 10'^'^ g and its main composition is 
hydrogen (70%), helium (29%) and less than 1% of more 
heavy elements, hke carbon, oxygen, etc. 

What are the nuclear processes which originate the 
huge thermonuclear energy of the Sun, and that has last 
4.6 X 10® years (the assumed age of the Sun)? It cannot 
be the simple fusion of two protons, or of a-particles, or 
even the fusion of protons with a-particles, since neither 
|He, |Be, or gLi, are stable. The only possibility is the 
proton-proton fusion in the form 



p + p 



d + c+ + Z/e 



(73) 



which occurs via /3-decay, i.e., due to the weak- 
interaction. The cross section for this reaction for pro- 
tons of energy around 1 MeV is very small, of the order of 
10"^'^ b. The average lifetime of protons in the Sun due 
to the transformation to deuterons by means of eq. ( 73 ) 



is about 10^° y. This explains why the energy radiated 
from the Sun is approximately constant in time, and not 
an explosive process. 

The deuteron produced in the above reaction is con- 
sumed almost immediately in the process 



P 



He 



-7- 



The resulting ^He reacts by means of 



me- 



He 



^Hc + 2p, 



(74) 



(75) 



which produces the stable nucleus |He with a great en- 
ergy gain, or by means of the reaction 



^He-I- ^He — > jBe -j- 7. 



(76) 



In the second case, a chain reaction follows as 



iBe 



2 (^He) 



(77) 



iBe 



|B +7, 2{iRc)+c++iy,. 

(78) 

The chain reaction (73)-(78) is called the hydrogen cy- 



cle. The result of this cycle is the transformation of four 
protons into an a-particle, with an energy gain of 26.7 
MeV, about 20% of which are carried away by the neu- 



trinos (see fig. 17 1. 

If the star contains heavier elements, another cycle can 
occur; the carbon cycle, or CNO cycle [37] • In this cycle 
the carbon, oxygen, and nitrogen nuclei are catalyzers of 
nuclear processes, with the end product also in the form 
4p — !• |He. fig. 18 describes the CNO cycle. Due to 
the larger Coulomb repulsion between the carbon nuclei, 
it occurs at higher temperatures (larger relative energy 
between the participant nuclei), up to 1.4 x 10*" K. In the 
Sun the hydrogen cycle prevails. But, in stars with larger 
temperatures the CNO cycle is more effective. 

After the protons are transformed into helium at the 
center of a star like our Sun, the fusion reactions start to 
consume protons at the surface of the star. At this stage 
the star starts to become a red giant. The energy gen- 
erated by fusion increases the temperature and expands 
the surface of the star. The star luminosity increases. 
The red giant contracts again after the hydrogen fuel is 
burned. 

Other thermonuclear processes start. The first is the 
helium burning when the temperature reaches 10® K and 
the density becomes 10^ g.cm"'^. Helium burning starts 
with the triple capture reaction 



3 (|He) 



7.65 MeV, 



(79) 



followed by the formation of oxygen via the reaction 



12 



C 



^He 



16 



+ 7. 



(80) 



For a star with the Sun mass, helium burning occurs 
in about 10^ y. For a much heavier star the temperature 
can reach 10® K. The compression process followed by the 
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FIG. 19: Relative distribution of elements in our galaxy. 



For temperatures above 3 x 10^ K more photo-nuclear 
processes appear. These yield more nuclei to be burned 
and heavier nuclei are produced: 



7 + 
^He + 



28 
14 
28 
14 



Si 
O 



?2Mg- 



He, 



llS+ 7, etc. 



(83) 



Due to the large number of free neutrons, many (n,7)- 
reactions (radiative neutron capture) elements in the 
mass range A = 28, . . . , 57 are formed. This leads to 
a large abundance of elements in the iron mass region, 
which have the largest binding energy per nucleon. For 
elements heavier than iron the nuclear fusion processes 
do not generate energy. 

For A > 100 the distribution of nuclei cannot be ex- 
plained in terms of fusion reactions with charged parti- 
cles. They are formed by the successive capture of slow 
neutrons and of /3~-decay. The maxima of the element 
distribution in = 50, 82, 126 are due to the small cap- 
ture cross sections corresponding to the magic numbers. 
This yields a trash of isotopes at the observed element 
distribution. 



burning of heavier elements can lead to the formation of 
iron and nickel. After that the thermonuclear reactions 
are no more energetic and the star stops producing nu- 
clear energy. 



B. Synthesis of heavier elements 

In figure [19] we show the relative distribution of ele- 
ments in our galaxy. It has two distinct regions: in the 
region A < 100 it decreases with A approximately like 
an exponential, whereas for A > 100 it is approximately 
constant, except for the peaks in the region of the magic 
numbers Z = 50 and N = 50, 82, 126. 



The thermonuclear processes (73 1- (80) can explain the 



relative abundance of ^He. 



^fC and gO. The processes 



occurring after |lle burning mainly form isotopes of igNe, 
12^S ilSi. We can understand the small abundance 
of the elements Li, Be and B as due to the small velocity 



with which they are formed via the reaction ( 76 1 and the 



first equation of (77), while they are rapidly consumed 



by the second reaction in (77) and the first reaction in 
(l78|. 



The synthesis of elements heavier than oxygen occur 
when, after helium burning, a new compression and heat- 
ing of the star rises the temperature to values higher than 
6 X 10® K. This situation ignites an intense carbon burn- 
ing: 



C 



c 



?°Ne- 



^He^ fJNa+ p 
n ItMg + 7. (81) 

Carbon and oxygen can also burn simultaneously: 



12c + ?tMg + |He, etc. 



16 



0+ 



?lSi- 



He, etc. 



(82) 



C. Thermonuclear cross sections 

The nuclear cross section for a reaction between target 
j and projectile k is defined by 



number of reactions target ^sec ^ ^/"j 
flux of incoming projectiles nkv 



(84) 



where the target number density is given by rij, the pro- 
jectile number density is given by Uk , and v is the relative 
velocity between target and projectile nuclei. Then r, the 
number of reactions per cm^ and sec, can be expressed 
as r = (Tvrijnk, or, more generally. 



-Vk\d rijdf 



rik- 



(85) 



The evaluation of this integral depends on the type of 
particles and distributions which are involved. For nuclei 
i and k in an astrophysical plasma, obeying a Maxwell- 
Boltzmann distribution (MB), 



2'KkT 



)3/2exp(- 



2kT ' ■ 



J' 



(86) 



eq. (85) simplifies to rj^k =< o'v > rijUk, where < av > 



is the average of av over the temperature distribution in 
p). 



In astrophysical plasmas with high densities and/or 
low temperatures, effects of electron screening become 
highly important. This means that the reacting nuclei, 
due to the background of electrons and nuclei, feel a dif- 
ferent Coulomb repulsion than in the case of bare nu- 
clei. Under most conditions (with non-vanishing tem- 
peratures) the generalized reaction rate integral can be 
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separated into the traditional expression without screen- 
ing and a screening factor 



IX. REACTIONS WITH RADIOACTIVE 
NUCLEAR BEAMS 



< (jv >*,fc= fscriZj,Zk, p, T, Yi) < av >j^k ■ (87) 

This screening factor is dependent on the charge of 
the involved particles, the density, temperature, and the 
composition of the plasma. Here Yi denotes the abun- 
dance of nucleus i defined by Yi — ni/(pNj^), where n,j 
is the number density of nuclei per unit volume and Na 
Avogadro's number. At high densities and low temper- 
atures screening factors can enhance reactions by many 
orders of magnitude and lead to pycnonuclear ignition. 

When in eq. (85) particle A: is a photon, the relative 



velocity is always c and quantities in the integral are not 
dependent on (Prij. Thus it simplifies to rj = Xj^-yUj and 
Aj.-y results from an integration of the photodisintegra- 
tion cross section over a Planck distribution for photons 
of temperature T 



7r2(c/i)" 



ca{E^)E': 
exp{KJkf) 



There is, however, no direct need to evaluate photodis- 
integration cross sections, because, due to detailed bal- 
ance, they can be expressed by the capture cross sections 
for the inverse reaction Z -I- m — >■ j -I- 7. 

A procedure similar to eq. (88) is used for electron 
captures by nuclei. Because the electron is about 2000 
times less massive than a nucleon, the velocity of the 
nucleus j is negligible in the center of mass system in 
comparison to the electron velocity (|uj — Ue| ~ I'^^eD- 
The electron capture cross section has to be integrated 
over a Boltzmann, partially degenerate, or Fermi distri- 
bution of electrons, depending on the astrophysical con- 
ditions. The electron capture rates are a function of T 
and rif, = Y^pN^, the electron number density [3^. In 
a neutral, completely ionized plasma, the electron abun- 
dance is equal to the total proton abundance in nuclei 

= ZiYi and rj= Aj,e(T, pYjn^. 

This treatment can be generalized for the capture of 
positrons, which are in a thermal equilibrium with pho- 
tons, electrons, and nuclei. At high densities {p > 10^^ 
g.cm~^) the size of the neutrino scattering cross section 
on nuclei and electrons ensures that enough scattering 
events occur to thermalize a neutrino distribution. Then 
also the inverse process to electron capture (neutrino cap- 
ture) can occur and the neutrino capture rate can be 
expressed similarly to Eq. ( 88 ) , integrating over the neu- 
trino distribution. Also inelastic neutrino scattering on 
nuclei can be expressed in this form. Finally, for nor- 
mal decays, like beta or alpha decays with half- life T1/2, 



we obtain an equation similar to Eq. ( 88 ) or rj of the 



last paragraph with a decay constant Xj = In 2/ Ti 



/2 



and 



The basic research activity in nuclear physics, driven 
by the desire to understand the forces which dictate the 
properties of nuclei, has spawned a large number of ben- 
eficial applications. Amongst its many progeny we can 
count reactor- and spallation-based neutron sources, syn- 
chrotron radiation sources, particle physics, materials 
modification by implantation, carbon dating and much 
more. It is an excellent example of the return to society 
of investment in basic research. 

All of these achievements have been realized by accel- 
erating the 283 stable or long-lived nuclear species we 
find here on Earth. In recent years, however, it has be- 
come evident that it is now technically possible to create 
and accelerate unstable nuclei and, and there are some 
6-7,000 distinct nuclear species which live long enough 
to be candidates for acceleration. They are the nuclei 
within the so-called drip-lines, the point where the nu- 
cleus can no longer hold another particle. This has led to 
many new opportunities in industry, medicine, material 
studies and the environment. 

Assume that a highly energetic uranium projectile 
{N/Z ^ 1.6) hits a target nucleus in an almost central 
collision, as shown in fig. [20j A part of the projectile 
(participant) is scrapped off and forms a highly excited 
mixture of nucleons with a part of the target. A piece of 
the projectile (spectator) flies away with nearly the same 
velocity of the beam. The neutron-to-proton ratio of the 
spectator is nearly equal to that of the projectile. Since 
the N/Z - ratio of light nuclei (stable) is close to one, 
the fragment is far from the stability line. Statistically, 
a large number of fragments with different N/Z - ra- 
tios is created and several new exotic nuclei have been 
discovered in this way. These nuclei can be collected in 
a secondary beam, further accelerated and induce reac- 
tions with a target nucleus. This method has become 
an important tool to study the properties of short-lived 
isotopes. 

Experiments with secondary-beams are limited by re- 
action cross section and luminosity. The luminosity L 
is defined as the product of beam intensity i and target 
thickness t, L = i ■ t. The reaction rate N is the product 
of luminosity and reaction cross section a^, N = ■ L.. 
In most of the reactions the usable target thickness is lim- 
ited by the width of the excitation function (i.e. the cross 
section as a function of the excitation energy). Produc- 
tion reactions with a wide excitation function covering a 
broad energy range can profit in luminosity by the use of 
thick targets. 

The condition for fragmentation of heavy ions is that 
the projectile should move faster than nucleons move in- 
side the nucleus. The projectile energy should be suffi- 
ciently above the Fermi domain, e.g. above 100 A MeV. 
The usable target thickness for these high energies is of 
the order of several grams per square centimeter, corre- 
sponding to 10^"^ atoms/cm^. The excitation function for 
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FIG. 20: (a) Schematic description of a nuclear fragmenta- 
tion reaction producing rare isotopes. The lower fragments 
are called participants, while the upper one is called the spec- 
tator. Using uranium projectiles {N /Z ~ 1.6) one expects to 
produce (light) spectator nuclei of about the same N/Z ra- 
tio, (b) Coulomb fission of relativistic projectiles leading to 
the production of rare isotopes. For a heavy unstable projec- 
tile an exchanged photon with the target can give it enough 
energy to fragment into several types of isotopes. 
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FIG. 21: Production cross sections for the tin isotopes from 
complete fusion (dotted line), fragmentation (solid line), and 
projectile fission of^^^V (dashed line). The symbols represent 
experimental data. The fragmentation cross-sections (solid 
line) have been calculated with a semi-empirical code [SueOO]. 



complete fusion of heavy ions, however, has a width of 
only 10 MeV. This corresponds to a usable target thick- 
ness of the order of one milligram per square centimeter 
or 10^^ atoms/cm^. Consequently beam intensities for 
the investigation of complete fusion reactions must be by 
four to five orders higher to achieve the same luminosity 
as for fragmentation. 

Fig. |2l] shows as an example the production cross sec- 
tions for the tin isotopes from complete fusion (dotted 
line), nuclear fragmentation (solid line), and Coulomb 
fission of ^ssy (dashed line). The symbols represent ex- 
perimental data. The fragmentation cross-sections (solid 
line) have been calculated with a semi-empirical code 
[SueOO]. 

The first experiments with unstable nuclear beams 
were designed to measure the nuclear sizes, namely the 
matter distribution of protons and neutrons. For stable 
nuclei such experiments are best accomplished with elec- 
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FIG. 22: The matter density radii of several light nuclei com- 
pared to the trend R ~ 1.18 A^^^ fm (dashed line) for 
normal nuclei. The solid lines are guides to the eyes. 



tron beams, which probe the nuclear charge (proton) dis- 
tribution. Electron scattering experiments with unstable 
beams can only be performed in an electron-nucleus col- 
lider. But the easiest solution is to measure the reaction 
cross section in collisions of unstable beams with a fixed 
target nucleus 38,. 

The reaction cross section in high energy collisions is 
given by 



aR = 2nJ[l- T(b)]bdb , 



(89) 



where 



Tib) = cxp 













— oo J 



PP (r)pT (R.-1-r) (fr 



(90) 

with R = (b, z) . (Tatjv is the nucleon-nucleon cross 
section at the corresponding bombarding energy, and 
Pp{T) is the projectile (target) matter density distri- 
bution. T(b) is known as the transparency func- 
tion. It is the probability that a reaction occurs for 
a collision with impact parameter b. The exponent 
[ctatat / pp (r )pt (R- -I- r) d^r] is interpreted as the 
mean-free path for a nucleon-nucleon collision. The re- 
action cross sections can be calculated using the matter 
distributions of the target in eq. [16] 

Experimentally it was observed a great increase in the 
rms radii f or t he neutron-rich isotopes ^He, ®He and ^^Li 
(see figure 22 1. Thus, the addition of the neutrons to ''He 
and ^Li nuclei increase their radii considerably. This can 
be understood in terms of the binding energy of the outer 
nucleons. The large matter radii of these nuclei have 
lead the experimentalists to call them "/la/o nuclei''^ . The 
binding energy of the last two neutrons in ^^Li is about 
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300 keV. In "^He it is 0.97 MeV. These are very small 
values and should be compared with S'„ =6 — 8 MeV 
which is the average binding of nucleons in stable nuclei. 
Abnormally large radii were also found for other light 
neutron-rich nuclei. 

The matter density radii of these nuclei do not follow 
the observed trend R ^ 1.18 A^^^ fm of normal nuclei. 
Thus the halo seems to be a common feature of loosely- 
bound neutron-rich nuclei. 

Several other methods have been devised to probe the 
structure of nuclei far from the valley of stability. Among 
these are (a) Coulomb dissociation [35^ , (b) Trojan horse 
method |40j . (c) asymptotic normalization coefficients 
[4T] . (d) heavy ion charge-exchange (42|, (e) knockout 
and (f) fusion reactions [44 . These methods yield 
different insights into the structure of exotic nuclei and 
comprise most of experiments in radioactive beam facil- 
ities l45l. 

Reactions producing rare nuclear isotopes has opened 
a new research front in nuclear physics with applications 
in many areas of science: (a) the rapid production of 



short-lived nuclei, many of which are poorly known, is 
important for astrophysics (r-process) and cosmology, (b) 
nuclear medicine benefits from studies of new nuclear iso- 
topes, (c) and the list goes on. But the basic question 
still remains: what combinations of neutrons and protons 
can make up a nucleus? The experimental detection of 
new nuclear isotopes is an ongoing research which will ul- 
timately lead to new insight and development of nuclear 
science, with enormous profit for mankind pUj . 
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